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We combine the ideas of scaling theory and universal conductance fluctuations with density- 
functional theory to analyze the conductance properties of doped silicon nanowires. Specifically, 
we study the cross-over from ballistic to diffusive transport in boron (B) or phosphorus (P) doped 
Si-nanowires by computing the mean free path, sample averaged conductance (G), and sample-to- 
sample variations std(G) as a function of energy, doping density, wire length, and the radial dopant 
profile. Our main findings are: (i) the main trends can be predicted quantitatively based on the 
scattering properties of single dopants; (ii) the sample-to-sample fluctuations depend on energy but 
not on doping density, thereby displaying a degree of universality, and (iii) in the diffusive regime the 
analytical predictions of the DMPK theory are in good agreement with our ab initio calculations. 

PACS numbers: 73.63.-b, 72.10.-d, 72.10.Fk 



Silicon nanowires (SiNWs) arc strong candidates for 
future nanoclcctronic and sensor applications 
most of the demonstrated devices the SiNWs are either 
p— or n— doped during the fabrication process. Very thin 
SiNWs with diameters below 5 nm have been synthesized 
by several groups @, @, @| , and due to the small cross sec- 
tion, scattering by dopants is likely to be very important. 
Moreover, due to reduced acoustic phonon scattering in 
quasi one-dimensional systems, long coherence lengths 
might be possible even at room temperature 0, thus 
emphasizing the importance of defect scattering. At the 
same time samplc-to-samplc variations become a crucial 
issue: when the device length and the mean free path 
are comparable, and shorter than the coherence length, 
variations of the positions of the individual dopant atoms 
can affect the conductance of the wire significantly. 

The mathematical theory of the conductance of disor- 
dered quasi one-dimensional systems has reached a high 
level of understanding in the diffusive as well as the local- 
ization regime @, Q. A standard way to model disorder, 
analytically as well as numerically, is to introduce ran- 
dom noise (Anderson disorder), which in a tight-binding 
description is included through randomly varying on-site 
energies. It does not seem obvious that a random disor- 
der is an adequate description of real physical disorder, 
such as dopants or vacancies, nor is there any obvious 
connection between a physical defect density and the am- 
plitude of the random disorder. In particular, if only a 
few dopants are present in a wire, the system is in the 
crossover region from ballistic to diffusive transport and 
the discrete and local nature of the impurities must be 
modeled adequately. 

Several recent theoretical works considered dopants in 
SiNWs using density-functional theory (DFT) [h], 
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12j , mainly focusing on the structural and energetic prop- 
erties of different radial dopant positions. As an impor- 



tant first step towards the modeling of physical SiNWs, 
Fernandez-Serra et al. [l3[ considered the scattering 
properties of single P dopants in thin nanowires. 

In this Letter we complete the analysis by calculating 
the conductance of long nanowires with a random dis- 
tribution of dopants (either P or B) along the wire. We 
calculate the sample averaged conductance (G). clastic 
mean free path (MFP) l e , localization length £, and the 
samplc-to-samplc fluctuations characterized by the sam- 
ple standard deviation std(G). We show that all these 
quantities can be understood and accurately estimated 
from the scattering properties of the single dopants, im- 
plying that relatively simple calculations arc sufficient in 
practical device modeling. The sample-to-sample fluc- 
tuations at a given energy and dopant type vary with 
L/l e in a universal way independently on the dopant con- 
centration, and in the diffusive regime for wire length 
l e < L < £, we observe good agreement with analyt- 
ical predictions of the Dorokhov-Mcllo-Pcreyra-Kumar 
(DMPK) theory [l5]. 

Method: The length and energy dependent conduc- 
tance of each sample is found using the Landaucr formal- 
ism together with a standard recursive Green's function 
(GF) approach where the full scattering region contain- 
ing the dopant atoms is constructed by repeatedly adding 
small unit cells [l6[ . The unit cells are constructed using 
first-principles local orbital DFT calculations 0, EH- 
We emphasize that our combined DFT and GF approach 
to calculate the conductance is fully ab initio within the 
low bias coherent transport regime. 

Single dopants: Figure [1] shows the transmission vs. 
energy through infinite hydrogen passivated wires con- 
taining only a single dopant atom placed at five different 
substitutional positions (1-5) indicated in the inset. The 
left part shows the transmission for B dopants at energies 
in the valence band while the right part shows results for 
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FIG. 1: (color online). Transmissions through wires with a 
single B (left) or P (right) dopant atom placed at the different 
radial positions shown in the inset. Left: energies relative to 
the valence band maximum, E v \ Right; energies relative to 
the conduction band minimum, E c . The band gap is 2.84 eV. 



sistance increases linearly for wire lengths shorter than 
the localization length, as shown in the left inset in 
Fig. [2] The initial linearly increasing resistance defines 
the energy dependent MFP, l e (E), through the relation 
R(L, E) = R c (E) + R c (E)L/l e (E) for wire lengths L < £, 
where R C (E) = 1/G (E) = h/(2e 2 N) is the contact re- 
sistance of a pristine wire with N conducting channels. 
It has recently been shown that this definition of the 
MFP agrees with values found using the Kubo formula 
and Fermi's golden rule HI, 21 1. In the linear resistance 
region, we suggest that on average, the scattering re- 
sistances from the dopants add classically according to 
Ohm's law, i.e. the mean resistance of a wire of length 
L < £ and average dopant-dopant separation d is given 
by 



(R(L, E)) = R C {E) + (R S (E)) L/d. 



(1) 



P dopants at energies in the conduction band. We no- 
tice that there are special resonant energies where the 
conductance abruptly drops, similar to those observed 
in Ref. due to enhanced local densities of states 

at the dopant atoms associated with quasi-bound states. 
As also pointed out in Ref. there are significant 

differences in the scattering properties between dopants 
located in the bulk of the wire and those situated at the 
surface. Interestingly, there are qualitative differences 
between P and B dopants. While B in position 1 (middle 
of wire) is a weak scatterer in the valence band, P at the 
same position is the strongest scatterer in the conduction 
band. For thicker wires the majority of the dopants are 
likely to be bulk-like and our calculations thus predict 
P-doped wires to have smaller mobilities than B-doped 
wires with the same dopant concentration, in agreement 
with experimental results (lij . 

Long wires: While the scattering properties of the 
single dopants are interesting in their own right, real 
nanowires contain several dopants with a certain distri- 
bution along the wire direction as well as in the radial 
direction. Due to interference effects between successive 
scattering events, it is not obvious if the single-dopant 
results carry over to the long wire case. In the rest of 
this Letter we investigate the conductance of long wires. 
In particular we examine to what extent the long wire 
properties, such as mean conductance (G) and variations 
std(G), can be determined from the scattering properties 
of the single dopants. 

A single recursive GF calculation yields the conduc- 
tance of a SiNW with a given dopant distribution, length 
and energy. At each energy we do this calculation for 300 
different realizations of the dopant positions, and we re- 
peat these calculations for a range of wire lengths, 10 
nm < L < 200 nm, all with the same dopant density. 
The dopants are distributed randomly along the wire di- 
rection [2^] as well as radially. The sample-averaged re- 



(R S (E)) is the average scattering resistance of the dif- 
ferent dopant positions which can be estimated from 
the single dopant conductances, G(E), in Fig. Q] as 
(R S (E)) = 1/(G(E)) - 1/G Q (E). From Eq. © we get 
an estimate of the MFP, l' e = R c d/(R S ). 
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FIG. 2: (color online). Energy dependent average conduc- 
tance calculated at L = 10 nm (circles) and L = 100 nm (dots) 
with the GF method. Solid lines follow from Eq. Q- Left 
inset: Resistance vs. wire length at energies E — E c = 0.01, 
0.06 and 0.16 eV (indicated with arrows in the main frame; 
same behavior is seen in the entire energy range). Solid lines 
are obtained using Eq. (TTJ. Right inset: (InT) vs. length at 
the same energies, where T — G/(2e 2 /h) is the transmission. 
Solid lines are linear fits in the interval 150 nm < L < 200 nm 
from which the localization length is determined. 

Figure [5] shows the sample averaged energy depen- 
dent conductance at wire lengths L = 10 nm (circles) 
and L = 100 nm (dots) for B and P doped wires with 
an average dopant-dopant separation d = 10 nm (corre- 
sponding to a bulk doping density of n w 10 19 cm~ 3 ) 
and uniform radial dopant distribution. The dashed line 
shows the conductance of the pristine wire, while the 
solid curves are estimated conductances obtained from 
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the single dopant transmissions shown in Fig. [T] using 
Eq. {T]). It is evident that the average conductances are 
well reproduced by the single-dopant results. 

We emphasize that Eq. (J]) is only valid in the quasi- 
ballisticc (L < l e ) and diffusive (l e < L < £) regimes. In 
the localization regime (L > £), the resistance increases 
exponentially, R(L) ~ exp(£/£). We have calculated the 
localization length from the (hiT) vs. L curves shown in 
the right inset of Fig. [2] as £ = — 1/slope, and the slope is 
determined from a linear fit in the interval 150 nm< L < 
200 nm. In Table U we show the resulting MFPs (l e ) and 
localization lengths (£) at four different energies. These 
are compared with the estimates l' e = R c d/(R S ) based 
on the single dopant transmissions, and £' = l' e (N + l)/2 
(this relation follows from random matrix theory Q). It 
is evident that both the MFP and the localization length 
can be estimated fairly accurately from the single dopant 
transmissions with a maximum error of 24%. We also 
note that the ratio t;/l e agrees with the prediction (TV + 
l)/2 with a maximum deviation of 30%. 
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TABLE I: MFP and localization lengths obtained from sample 
averaging (l e and £), and estimated values, l' e and obtained 
from the single dopant transmissions. TV is the number of 
conducting channels. 



Sample fluctuations: An interesting question is 
whether or not the sample-to-sample variations can also 
be understood in terms of the single-dopant transmis- 
sions. In the diffusive regime we expect the variations to 
be close to the universal conductance fluctuation (UCF) 
value, 0.73 e 2 /h for quasi-one dimensional systems [22I ]. 
Figure [3] shows the standard deviation, std(G), plotted 
against normalized length, L/l e , for B (a)-(c) and P 
(d)-(f) doped wires. The conductance fluctuations cor- 
responding to different concentrations lie close to each 
other, and we therefore conclude that at a specific en- 
ergy and type of dopant, the sample variations are in- 
dependent of dopant concentration but only depend on 
the normalized length. This is in accordance with the 
theory of UCF [22| and single parameter scaling theory 
|23j ] which predict that the sample fluctuations are inde- 
pendent on the disorder strength. For modeling purposes 
this is a very convenient result because one can limit the 
simulations to only one dopant concentration. 

In the diffusive limit, L/l e > 1, analytical results are 
known for std(G(£)). For quasi one-dimensional systems 
with many conducting channels and weak random dis- 
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FIG. 3: (Color online) Sample standard deviation, std(G) 
vs. normalized length L/l e with average dopant-dopant 
separation d = 19 nm (dots), d = 10 nm (squares) and 
d — 5nm (circles). (a)-(c) correspond to the energies 
E — E v = -0.12, -0.20, -0.30 eV in the valence band for 
B doped wires, (d)-(f) correspond to the energies E — E c = 
0.06, 0.16, 0.26 eV in the conduction band for P doped wires. 
The solid lines are analytical solutions to the DMPK equation 
of Ref.0 which equals the UCF value 0.73 e 2 /h at L = 0. 



The small arrows mark the maximum std(G), 
from the single-dopant transmissions. 



estimated 



order, the DMPK equation [l 
dependence, std(G(L)) 



predicts a weak length 
>/8/15 - 64/315 • £/£ • e 2 /h, 
shown in Fig. [3] (solid lines). As the length of 
the wire decreases, the sample fluctuations increase, and 
in most cases a maximum in std(G) is reached around 
L/l e = 0.5. In the limit L — > there are no dopants in 
the wire and the sample-to-sample variations vanish. 

We next address this maximum value and the general 
behavior of the sample-to-samplc fluctuations vs. L using 
the single-dopant conductance results. Figure SJa) shows 
the maximum values of std(G) vs. energy for P-doped 
wires with a uniform radial distribution (squares) and 
pure surface doping (circles). In the former case all the 
dopant positions 1-5 are equally probable[24[ while in 
the latter the dopants can only be in position 4 or 5, cf. 
inset in Fig. [TJ The lower solid line in FigQJa) shows 
the standard deviation s =std({Go, Go, G4, G5}), where 
Go is the conductance of a pristine wire while G4 and G5 
are the single-dopant conductances with the P dopant 
placed at position 4 and 5. This sequence represents a 
situation where there is 50% chance for a pristine wire, 
and 50% chance for a wire with a single dopant either at 
position 4 or 5. The upper solid line in Fig. |4ja) shows 
similar values for the uniform dopant distribution. When 
s is larger than the UCF value (the horizontal dashed 
line), the maximum std(G) clearly follows the trends in 
s. When s is smaller than the UCF value, the maximum 
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std(G) lies close to the UCF level. 
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FIG. 4: (color online). Maximal standard deviation (a) for 
for P doped wires (d — 10 nm,) with a homogeneous radial 
distribution (squares) and a pure surface doped wire (circles). 
The solid lines show the standard deviation, s, among the 
single dopant transmissions and the pristine wire, and the 
dashed line marks the UCF value 0.73 e 2 /h. Panel (b) shows 
the MFP vs. energy for the two radial distributions (circles 
and squares). The lines represent values obtained from the 
single dopant transmissions. 

We conclude that the shape of the std(G) vs. L/l e 
curves in Fig. [3J can be qualitatively predicted from the 
single-dopant transmissions: When s, indicated with 
small arrows in Fig. [3J is smaller than the UCF value 
(0.73e 2 //i), we expect std(G) to approach the analytical 
line from below. The B doped wires in Fig. [3jb) show 
such a behaviour. Otherwise, if s is larger than the UCF 
value, std(G) will reach a maximum value close to s and 
approach the analytical line from above. 

Mean free path: Figure[3Ja) shows that the sample fluc- 
tuations are significantly reduced in the surface doped 
wires due to the weaker and more uniform scattering 
properties of the surface positions. The MFP also de- 
pends significantly on the radial distribution as seen in 
Fig. EJb), showing the MFP for P doped wires with a 
uniform (circles) and pure surface (squares) distribution 
of dopants. The solid lines are estimated from the single- 
dopant transmissions using Eq. ([1]) once again showing 
that the mean values of the conductance can be accu- 
rately estimated from the single dopants. A very signif- 
icant increase in the MFP is observed for the surface 
doped wires as compared to the uniform distribution. 
This might suggest that increased device performance 
could be achieved if the P dopants are located close to 
the surface, which indeed are the energetically most fa- 
vorable positions in the wires studied here and also found 
m Ref. Hiim. There are, however, potential problems 
with dopants close to the surface, as they can be passi- 
vated by addition of an extra hydrogen atom [l2|, EH ■ 

In conclusion, we have considered the conductance 



properties in B and P doped SiNWs. We find that the 
sample averaged conductance and the samplc-to-samplc 
fluctuations as well as the mean free path and localiza- 
tion length can be predicted quantitatively from the scat- 
tering properties of the single dopants. Time consuming 
sample averaging can thus be avoided, which greatly sim- 
plify modelling of the statistical conductance properties. 
These findings may have a high impact on first-principles 
modelling of electron transport in nanowires. 
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